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ABSTRACT 



We determine the ratio of helium- to hydrogen-atmosphere white dwarf stars 
as a function of effective temperature from a model atmosphere analysis of the 
infrared photometric data from the Two Micron All Sky Survey combined with 
available visual magnitudes. Our study surpasses any previous analysis of this 
kind both in terms of the accuracy of the T e g determinations as well as the size 
of the sample. We observe that the ratio of helium- to hydrogen-atmosphere 
white dwarfs increases gradually from a constant value of ~ 0.25 between T e g = 
15, 000 K and 10,000 K to a value twice as large in the range 10, 000 > T cS > 8000 
K, suggesting that convective mixing, which occurs when the bottom of the 
hydrogen convection zone reaches the underlying convective helium envelope, is 
responsible for this gradual transition. The comparison of our results with an 
approximate model used to describe the outcome of this convective mixing process 
implies hydrogen mass layers in the range Mn/M to t = 10~ 10 to 10~ 8 for about 
15% of the DA stars that survived the DA to DB transition near T e g ~ 30, 000 K, 
the remainder having presumably more massive layers above Mft/M tot ~ 10~ 6 . 

Subject headings: Convection — stars: atmospheres — stars: interiors — white 
dwarfs 



1. INTRODUCTION 

The Two Micron All Sky Survey (2MASS) represents one the largest set of homogeneous 
infrared photometric data for all areas of the sky applicable to white dwarf stars of various 
spectral types. Indeed, the 2MASS survey contains JHKs magnitudes for almost all point 
sources of the sky up to the magnitude thresholds of J ~ 16.8, H ~ 16.5, and Ks ~ 16.2 for 
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detections with any formal uncertainty (see for instance Fig. 1 of lTremblay fc Bergeronll2007l ). 
This yields a nearly complete magnitude limited survey of white dwarfs that allows us to trace 
a portrait of the local population of white dwarfs, assuming we can properly identify the white 
dwarfs in this survey. To that effect, the most complete compilation is that of the Villanova 
White Dwarf Catalog^ (WDC), which contains more than 5500 spectroscopically identified 
white dwarfs reported in the literature. The interest of such a sizeable sample of homogeneous 
infrared photometric observation s is to reach an unprecedented level in the statistical analysis 
of white dwarf stars. In addition, iTremblay fc Bergeron! (120071 ) have demonstrated the overall 
reliability of the 2MASS Point-Source Catalog (PSC) photometry for the purpose of detailed 
comparisons with the predictions from white dwarf model atmospheres. 

Studies of the spectral evolution of white dwarfs attempt to understand how the surface 
composition and the corresponding spectral s i gnatu re evolve along the white dwarf cool- 
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1991 for a review). Because of the intense 



mg sequence (see 

gravitational field present in these stars, hydrogen will float to the surface in the absence 
of competing mechanisms, while heavier elements will sink below the photospheric regions. 
Indeed, the majority of white dwarfs have optical spectra that are completely dominated by 
strong hydrogen Balmer lines — the DA stars — with hydrogen-dominated atmospheres. 
We also find that almost 25% of white dwarfs have helium-dominated atmospheres, mostly 
DO, DB, DQ, and DZ stars, and some of the DC stars as well. The fact that not all white 
dwarfs have hydrogen-rich atmospheres imply that either some stars are born with hydrogen- 
deficient atmospheres and remain that way throughout their evolution, or alternatively, that 
several physical mechanisms such as diffusion, accretion, convection, radiation pressure, and 
stellar winds are competing with gravitational settling in determining their surface compo- 
sition as they cool off. This particular issue can be resolved by analyzing the distribution of 
spectral types as a function of effective temperature (or any tem perature in dex) for a large 
sample of white dwarf stars, in a way similar to the analysis of ISionl (11984 ) nearly twenty 
five years ago. 

Before discussing the physical processes that will affect the spectral evolution of hydrogen- 
and helium-atmosphere white dwarfs, it is important first to consider the constraints brought 
about by the studies of the hottest white dwarfs. The canonical mass fractions of light el- 
ements in white dwarfs considering the mass threshold for residual nuclear burning are 
M He /M tot ~ 10~ 2 and M H /M tot ~ 10" 4 . However, multiple excursions on the AGB nu- 
clear burning phase, during which ve ry late helium flashes could remove essentially all of 
the hydrogen (jWerner fc Herwigj|2006l ). could also produce white dwarfs with total hydrogen 
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masses much lower than the canonical value of 10~ 4 . The first class would correspond to the 
progenitors of the hottest DA stars while the second class would represent the progenitors 
of hot white dwarf stars with hydrogen-deficient atmospheres — the PG 1159 stars at very 
high effective temperatures, and the DO stars at slightly lower temperatures, whose spectra 
are dominated by Hell lines. It has traditionally been believed that by the time DO stars 
cool off to T e ff ~ 45,000 K, residual amounts of hydrogen (Mn/M to t ~ 10~ 16 ) thoroughly 
mixed in their helium-rich envelope would gradu ally accumulate to the surface, tur ning all 
helium-atmosphere white dwarfs into DA stars OFontaine Sz Wesemaell 119871 . Il99ll ). This 
scenario would account for the existence of the so-called DB gap, a range of temperature 
between T e fj = 45, 000 K and 30,000 K where all white dwarfs have a DA spectral type. 
However , the recent discovery i n the Sloan Digital Sky Survey of several DB stars within 
the gap (jEisenstein et al.ll2006al ) implies that the total mass of hydrogen left in the envelope 
of DO stars can be even smaller than previously believed. It is then generally admitted 
that the significant increase in the number of helium-atmosphere DB white dwarfs below 
T e fj ~ 30, 000 K can be explained in terms of the convective dilution of the superficial hydro- 
gen atmosphere by the underlying helium convective envelope, provide d that the hydrogen 



layer is sufficiently thin (M H /M tot ~ 10~ 15 . MacDonald fc Venneslll99lh . 



At even lower effective temperatures (T efT < 12,000 K), hydrogen-atmosphere white 
dwarfs get a second opportunity to evolve into stars with helium-dominated atmospheres 
when the superficial hydrogen layer becomes convective over a significant fraction of its 
depth. This process is illustrated in Figure [T] where the extent of the hydrogen convection 
zone is displayed as a function of decreasing effective temperature for a 0.6 M Q DA white 
dwarf, based on evolu tionary models with thick hydrogen layers similar to those described by 
Fontaine et al.l (120011 ) and kindly provided to us by G. Fontaine and P. Brassard. These cal- 
culations show that if the hydrogen envelope is thin enough, the bottom of the hydrogen con- 
vection zone may eventually reach the underlying and more massive convective helium layer, 
resul t ing in a mixing of the hydrogen and h e lium layers (IStrittmatter fc Wickramasinghe 
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7J). Fig- 



197ll : IShipmaiJll972lBaglin fc Vauclair! 1 19731 : iKoesteri Il97a 
ure [1] also indicates that the effective temperature at which this mixing occurs will depend 
on the thickness of the hydrogen envelope. The thicker the envelope, the lower the mixing 
temperature; if the hydrogen layer is more massive than Mn/M tot ~ 10 -6 , mixing will never 
occur. 

The simplest physical model that can be used to describe this convective mixing pro- 
cess is to assume that hydrogen and helium are homogeneously mixed. Since the helium 
convection zone is much more massive (MH e -conv/-^tat ~ 10~ 6 ) than the hydrogen layer 
when mixing occurs, it is generally assumed that a DA star would be transformed into a 
helium-atmosphere white dwarf of type DB, DQ, DZ, or DC with only a trace abundance 
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of hydrogen. More detailed calculations discussed by iFontaine fc Wesemaell (ll99ll ) confirm 
these predictions. If such a process takes place in cool white dwarfs, we then expect the ratio 
of helium- to hydrogen-atmosphere white dwarfs to increase at low effective temperatures. 
A comparison between the observed ratio as a function of T e g- and the theoretical expecta- 
tions would make it possible to estimate the thickness of the hydrogen layers in DA white 
dwarfs. In turn, these determinations could be compared to independent measurements of 
the hydrogen layer mass inferred from ZZ Ceti asteroseismology. 

The best approach to study the occurrence of this convective mixing process begins with 
the statistical analysis of large samples of cool white dwarfs for which we can determine the 
main atmospheric constituent and the effe ctive t empe rature, or any analogue temperature 
index such as the absolute V magnitude. ISionl (119841 ) was the first to study this problem 
by estimating the non-DA to DA ratio at low effective temperatures using a proper motion 
sample of 695 white dwarfs drawn from the catalog of s pectr o scopic ally identified white 
dwarfs available at that time. The results from Figure 1 of ISionl ( 19841 ) are reproduced here 
in Figure [2] in terms of the non-DA to DA ratio (rather than the absolute number of DA and 
non-DA stars) as a function of the absolute visual magnitude My. Looking at the results, we 
notice a first increase in this ratio at My > 11.25, which corresponds to T e g ~ 15, 000 K for a 
0.6 M Q white dwarf, a temperature much in excess of the theoretically predicted value for the 
convective mixing process. A second increase occurs at My ~ 12.5, or T e g ~ 10000 K, but 
the observed ratio drops suddenly afterwards and this increase is probably not significant. 
Then a third increase occurs at My > 13, or T cff < 8000 K. Clearly, the evidence based 
on these results are not exactly convincing, and our current picture of the situation below 
12,000 K is at best s ketchy. The global portrait did not change significantly since the 
analysis of ISionl (119841 ). and this remains the only available proof quoted in the lit erature of 
the transformation of some DA stars into non-DA stars (see also lGreensteinlll986l ). 



Our goal is to improve upon the analysis of ISionl (119841 ) by using the 2MASS photo- 
metric sample combined with the WDC database to identify white dwarfs and determine 
their atmospheric composition, and most importantly to obtain more accurate temperature 
determinations using a full model atmosphere analysis of the VJHKs photometry. We first 
describe in § 2 the 2MASS photometric white dwarf sample used in our analysis. The de- 
termination of the atmospheric parameters is then presented in § 3 and the uncertainties 
related to our approach are discussed at length in § 4. The evolution of the ratio of helium- to 
hydrogen-atmosphere white dwarfs as a function of effective temperature is then determined 
in 5 5. Our conclusions follow in § 6. 
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2. DEFINITION OF THE 2MASS WHITE DWARF SAMPLE 

Our starting point is the Villanova White Dwarf Catalog (WDC) of spectroscopically 
identified white dwarfs. We took the most recent electronic version of the catalog (August 
2006), which includes a total of 5557 white dwarfs. We made a first selection by taking only 
the objects with at least one published visual magnitud^l and with at least one measurement 
brighter than V = 18.5. Figure [3] shows the predicted V magnitude for which .J^mass = 16.5 
- a representative value for the 2MASS detection limit in this band — as a function of T e g 
for our pure hydrogen and pure helium model atmospheres at logg = 8. Only white dwarfs 
located below these curves at a given effective temperature can be detected by 2MASS in 
the J and H bandpasses. For white dwarfs with T e g- > 5000 K, the V < 18.5 criterion 
is sufficient to find a maximum number of white dwarfs. Finally, since our study is aimed 
at tracing the spectral evolution of cool white dwarfs, we neglected all objects with early 
spectral classifications (e.g., PG 1159, DO, etc.) or for which the spectral classification was 
unavailable altogether. Thus, a total of 1473 remaining white dwarf stars define our initial 
WDC sample. 

We then searched the 2MASS PSC for JHKs photometric observations of all white 
dwarfs from this initial WDC sample. We used the GATOR batch file tool and a 20" search 
window centered on a set of improved coordinates measured by J. B. Holberg (2005, private 
communication). In most instances, multiple sources were found within the search window 
and we unambiguously identified each object by comparing the 2MASS atlas with the finding 
charts available from the online version of the Villanova WDC. For several objects close to 
the detection limits, only one or two magnitudes are detected with formal uncertainties. We 
neglect all objects with only one detected magnitude since in such cases, only one degree of 
freedom would be available in our minimization procedure described in § 3. We recovered 
940 objects in the 2MASS PSC from the 1473 objects contained in our initial WDC sample. 
The remaining 533 objects are eliminated for the following reasons: 508 are too faint to be 
detected by 2MASS in at least two bands, 15 are bright enough in the 2MASS atlas but 
are not part of the 2MASS PSC, and 10 could not be unambiguously identified from the 
comparison of the 2MASS atlas and the published finding charts. 

A spectral energy distribution is then built for each object by combining the 2MASS 
JHKs photometry with the V magnitude (or Stromgren y) taken from the W DC. In addi- 



tion, about 10% of our V magnitudes come from unpublished USNO photometry fjDahn et al. 



20071 ). When more than one visual magnitude is available, we take an unweighted average, 



2 We have neglected magnitudes prior to 1972, photographic magnitudes, B magnitudes, as well as those 
flagged as uncertain. 
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unless we had reasons to believe that a given measurement was less accurate than others. 
We also assume a typical uncertainty of 0.05 mag for the visual magnitudes. 



ATMOSPHERIC PARAMETER DETERMINATIONS 



We classify the white dwarfs in our sample into two broad categories, those with hydro- 
gen atmospheres (DA stars, including DAB and DAZ stars) and those with helium atmo- 
spheres (DB, DZ, and DQ stars as well as helium-atmosphere objects with weak hydrogen 
features such as Ross 640 and L745-46A; the case of DC stars is discussed further below). 
When the spectral classification is ambiguous, we check in the literature to confirm the pres- 
ence of hydrogen lines compatible with the estimated effective temperature. We then fit 
the VJHKs (or VJH) p hotometry with our gr i ds of p ure hydrogen or pure helium model 
atmospheres described in iTremblav fe Bergeron! (120071 ) and refere nces therein. Our fitting 
technique is similar to that outlined in Bergeron et ali ( 1997 . 2001 ). Br iefly, the magnitudes 
are co nverted into monochromatic fluxes using the zero points defined in iHolberg fc Bergeron 
(120061 ) for photon counting devices, but using the transmis s ion fu nction of iBessell ( Il990h for 
V and the 2MASS transmission functions of ICohen et al.l (120031 ) for JHKs- The resulting 
fluxes are then compared with those predicted from the model atmospheres, properly av- 
eraged over the same bandpasses. The effective temperature and the solid angle ir(R/D) 2 
(with R the radius of the star and D the distance to Earth) are considered free parame- 
ters in the \ 2 minimization procedure. We assume a value of log g = 8.0 throughout. The 
surface gravity has very little effect on the predicted (V — J) color indices of white dwarfs. 
The results shown in Figure H] indicate that a difference of 1 dex in log g corresponds to a 
color difference smaller than 0.05 mag at all temperatures. This slight difference, which is 
of the order of our color uncertainties, allows us to fix the surface gravity at log g = 8.0 for 
all objects. Since the mass distributions of hydroge n- and helium- atmosphere white dwarfs 



have to a good approximation the s ame mean value (IBergeron et al.lll992 
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Liebert et al . 2005; Voss etal 1 l2007h . this simplifying assumption should not affect the results 
of our analysis significantly. 

For the featureless DC stars, we consider both hydrogen- and helium-atmosphere solu- 
tions. If the hydrogen solution is above T c g = 5000 K, we adopt the helium solution since 
in principl e, Ha should have been det ected spectroscopically according to the photometric 
analyses of IBergeron et al.l (119971 . l200ll ). However, if the hydrogen solution is below this tem- 
perature, it is impossible to determine the main atmospheric constituent from the restricted 
available photometry, and therefore we set the limit of our analysis to T e g > 5000 K. 



There is also a significant fraction of binary white dwarfs in our sample. If the com- 



- 7- 



panion is a cool main- sequence star, the observed JHKs fluxes are substantially larger than 
those predicted by the mod els, and the (J — H) color index takes a large positive value 



( ITremblay fc Bergeronll2007l ). which typically yields a bad fit to the energy distribution. We 
eliminated 90 binary candidates satisfying at least two of the following criteria^]: (1) a pho- 
tometric temperature below 6000 K incompatible with the spectral type, (2) an observed 
flux at H larger than that at J, and (3) a flux at J larger than ten times the flux at V. We 
have also eliminated 9 binary systems that are partially resolved in the 2MASS atlas, but 
whose colors in the the 2MASS PSC are likely to be contaminated. The large majority of 
our binary candidates have been previously identified as such in the literature. 

Sample fits for eight bright (11 < V < 13.5) and well studied white dwarfs in our 
sample are displayed in Figure [5J We notice in each case that the best fit is in excellent 
agreement with the V JHKs photometry, for both hydrogen- or helium-atmosphere white 
dwarfs. We also compare in Table 1 the photometric temperatures for these objects with 
independent spectroscopic or photometric effective temperature determinations reported in 
the literature. This comparison indicates that our technique applied to cool white dwarfs 
of various spectral types yields T e g values with a precision of the order of 5% compared to 
other methods, provided we have accurate photometric data. 



4. ATMOSPHERIC PARAMETER UNCERTAINTIES 

We attempt in this section to better quantify the uncertainties of our effective temper- 
ature determinations in order to determine the optimal bin widths for the histograms used 
in our study of convective mixing, presented in section § 5. These uncertainties can be first 
quantified by looking at the rate of change of the (V — J) color index as a function of T e g. 
We show in Figure [6] the effective temperature variation, AT e g, for a corresponding change 
of the (V — J) color index of A(V — J) = 0.08, a value that defines the mean uncertainty of 
this color index in our sample. We notice that the T e g- uncertainties increase rapidly towards 
high effective temperatures, and for this reason, we restrict our analysis to T c g < 15, 000 K. 

We compare in Figure [71 o ur photometric T R ff de t ermin ations for 133 cool white dwarfs in 



common with the analyses of iBergeron et al.l (119971 . |200l|) , which are based on independent 



BVRIJHK photometric data and in several cases, trigonometric parallax measurements 
from which surface gravities can be determined. To be internally consistent, we have revised 
their atmospheric parameters using our own fitting procedure based on the slightly different 



3 16 of these did not exactly meet two criteria but they have poor fits and atmospheric parameters that 
are incompatible with the published spectral type. 
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synthetic photometric calibration of iHolberg fc Bergeron! (120061 ). As well, we consider only 



the stars that have been observed in at least the J and H bands. We can see that the agree- 
ment is excellent in this temperature range and that the VJHKs photometry is sufficient 
to constrain the T e g values with the precision sought in our analysis. 

We present in Figure [8] a similar comparison of our photometric temperatures with in- 
dependent values obtained from line profile fitting of DA and DB stars. The spectroscopic 
temperatures for the 198 DA stars in common with our sample are taken from the ongo- 
ing spectroscopic survey of A. Gianninas (2007, private communication). We notice in this 
case that both samples are generally in good agreement, but several discrepancies are also 
observed where the spectroscopic temperatures are significantly larger than the photometric 
values (AT c ff > 0), in particular at high temperatures. A first explanation for this discrep- 
ancy is that above T e g- > 11,000 K, the (V — J) color index is a slowly varying function 
of effective temperature (see Fig. [6]), and thus white dwarfs with large temperature uncer- 
tainties contaminate this particular region. This is confirmed by the fact that most white 
dwarfs showing a large discrepancy are near the 2MASS detection limit (objects not detected 
at Ks in 2MASS are shown by open circles in Fig. [8]). Therefore, we add the constraint 
that only white dwarfs detected in all three JHKs bands are allowed in our analysis for 
T eff > 11, 000 K. The second group of objects showing a positive AT eff are unresolved double 
degenerate systems, three of which are identified in Figure [8] by filled squares. In all cases 
the photometric temperature is substantially lower than the spectroscopic value, and both 
spectroscopic and photometric fits are excellent. It is therefore nearly impossible to eliminate 
double white dwarf systems from our analysis. There remain 10 other objects with a AT c g 
larger than 20% for which the photometric fit is good and there is no clear explanation for 
a discrepancy. These could also correspond to unresolved degenerate binaries, although a 
closer inspection rather suggests that the V magnitudes, which are old measurements, are 
probably in error. 

In the bottom panel of Figure [8] we compare our photometric t e mperatures for 18 



DB st ars with those determined spectroscopically by iBeaucharnpl (119951 ); iBeauchamp et al. 
( jl996l ). We see that the agreement is excellent although the spectroscopic temperatures are 
systematically larger than the photometric values by about 5%. This small shift could be 
due to the neglect of small traces of hydrogen or metals in both analyses, or to inaccuracies 
in the treatment of the broadening theory for neutral helium lines at lo w effective tempera - 



tures, which will mainly affect the spectroscopic temperatures (see, e.g.. lKepler et al.ll2007l ). 
The effective temperature dispersions observed in the figures above indicate that histograms 
with bin widths of 1000 K below T g = 11, 000 K and of 2000 K above this temperature are 
acceptable. These values are reproduced in Figures [6] to [HJ 
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RESULTS ON CONVECTIVE MIXING 



We present in Table 2 the number of hydrogen- and helium-atmosphere white dwarfs 
in different effective temperature bins for our final 2MASS sample of 447 objects (after the 
temperature cutoff at T c g = 15, 000 and 5000 K and the removal of 2MASS data with low 
accuracy above T eS = 11,0 00 K have b een applied). This sample is nearly twice as large as 
that used in the analysis of ISionl (119841 ) in the same range of effective temperature. We also 
display in Figure [9] the corresponding histogram of the number of hydrogen- and helium- 
atmosphere white dwarfs as a function of T e g. The total number of stars in each bin below 
T e fj = 11,000 K is comparable, a result that can be explained by the competition between 
the increasing luminosity function and the decreasing detection probability of cooler, fainter 
objects. We further observe that hydrogen-atmosphere DA white dwarfs are always the 
dominant type. Consequently, convective mixing — if this process occurs at all — should 
affect only a small fraction of cool DA stars. 

For a magnitude limited sample — here the H magnitude — , we probe a limited volume 
of space. This volume is defined in terms of the distance at which a given white dwarf can 
be detected at the limiting magnitude of the survey. However, this distance is different 
for hydrogen- and helium-atmosphere white dwarfs because they have different continuum 
opacity sources and corresponding monochromatic fluxes. Figure [TU1 presents the ratio of the 
volume probed by hydrogen and helium white dwarfs as a function of effective temperature 
for a magnitude limit in the H bandpass. While the assumption of a fixed value of log g 
was appropriate for determining the effective temperatures, this may not be the case for the 
volume effect discussed here. Indeed, white dwarfs stars with large radii (or low masses) 
can be detected at larger distances. The differ ences in the mass distri butions of DA and DB 
stars could therefore aff ect the ratio (see also lEisenstein et al.ll2006al ). While the results of 
Beauchamp et al.l (119961 ) suggest t hat the mas s distri bution of DB stars is narrower than that 
of DA stars, the recent results of IVoss et al.l (120071 ) reveal instead that both distributions 
are comparable. Hence the small residual differences in the low and high mass tails are not 
expected to affect our results significantly. 

To take into account the volume effect, we henceforth correct the ratio of the number of 
helium- to hydrogen- atmosphere white dwarfs, Ane-atm/An-atm, by the corresponding factor 
given in Figure 11. This ratio corrected for the volume effect is reported in Table 2 and 
displayed in Figure [IT] as a function of effective temperature. The error bars are statistical 
and they are mainly governed by the limited number of helium-atmosphere objects in each 
bin. We also show in the same figure the ratio obtained by shifting the temperature bins in 
Figure! by 500 K. 



We can immediately notice one important trend in the observed ratio of helium- to 
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hydrogen-atmosphere white dwarfs displayed in Figure [TTJ At high effective temperatures, 
this ratio takes a constant value of ~ 0.25 and it shows a gradual increase from T e fj ~ 
10, 000 K to a fairly constant value twice as large (~ 0.5) at 8000 K. The only physical 
mechanism that could account for this increase is the convective mixing of thin hydrogen 
layer with the more massive underlying helium envelope. These results further imply that 
about 15% of all DA stara^ between T e g = 15, 000 K and 10,000 K have turned into helium- 
atmosphere white dwarfs — either DQ, DZ (or DZA), or DC — by the time they reach 
8000 K. Note that the results of our analysis become too uncertain above the temperature 
range displayed here (see Fig. [6]). 

We can gain a better understanding of our results by looking at the simple theoretical 
model of convective mixing discussed in the Introduction. As the white dwarf cools off, 
the bottom of the hydrogen convection zone will eventually reach the helium convection 
zone, provided the hydrogen layer is thin enough (see Fig. [I]), with the mixing temperature 
depending on the thickness of the hydrogen layer. White dwarfs with thicker hydrogen layers 
will mix at lower effective temperatures. A comparison of the results shown in Figures [1] 
and [TT] implies that there are very few DA stars between T c g ~ 15, 000 K and 10,000 K 
- the region where ZZ Ceti stars are found — with hydrogen layer masses smaller than 
Mn/Mtot ~ 10~ 10 . The gradual increase of the observed ratio from 10,000 K down to 8000 K 
further implies that only 15% of DA white dwarfs in the 10—15,000 K temperature range have 
superficial hydrogen layers with a mass in the range Mn/Mtot ~ 10 -10 to 10 -8 . Note that this 
range of layer mass depends sensitively on the assumed convective efficiency (ML2/a = 0.6 
in Fig. [T]). Envelope models with a more efficient parameterizatio n of the convective e nergy 
transport would imply even thicker hydrogen layers (see Fig. 9 of iTassoul et~aDll99oh . The 



lack of evidence for convective mixing at lower effective temperatures between T e g = 8000 K 
and 5000 K implies a discontinuity in the hydrogen layer mass distribution, with the large 
majority (~ 85%) of all DA stars between T eS = 15, 000 K and 10,000 K having hydrogen 
layer masses above Mn/Mtot ~ 10~ 6 or so. The constant ratio of ~ 0.25 in this temperature 
range also implies that ~ 20% of all DA stars in the DB gap (or DB "trough") between 
T e ff = 45, 000 K and 30,000 K have turned into DB stars by the time they reach 15,000 K; 
we cannot tell from our results at which effective temperature the DA to DB transition 
occurs precisely. 

The numbers given above could change slightly if one takes into account the presence 
of traces of metals (including carbon) and hydrogen in the model atmosphere calculations of 
DQ and DZ stars. Indeed, the atmospheric structure of cool helium-atmosphere white dwarfs 



4 The transformation implies the following change of ratio A^Ho-atm/^VH-atm — * (-^He-atm 

—transformed — atm ransformcd ) • 
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is affected signi ficantly by small traces of metals or hydrogen, w hich then provide most of the 
free electrons (IProvencal et al.l 120021 ; iDufour et al.l 120051 . 120071 ). These additional electrons 
increase the contribution of the He~ free-free opacity, which in turn reduces the atmospheric 
pressure. The results of Dufour et al. (2005, 2007) reveal that photometric temperatures 
of DQ and DZ stars can be systematically reduced by 500 K on average when metals are 
included in the model calculations. Therefore, this could affect the precise temperature 
dependence of the change in the ratio of helium- to hydrogen-atmosphere white dwarfs 
observed in Figure [TTJ 

We finally note that our results are compatible with the asteroseismological analyses 
of ZZ Ceti stars, which suggest that most ZZ Ceti st ars have fairly thick hydrogen layers 
of th e order o f Mjj/M tnt ~ 10~ 6 and higher (see, e.g.. iFontaine et al.l 19921 ; iBergeron et al. 



19931 : lBradlevlll998l : Brassard fc Fontaine! 12009 : IPech fc Vauclair! 120061 ) . Only a small frac- 
tion (about 15%) of all ZZ Ceti stars should have thin layers according to our results. 
Furthermore, our re sults are also in qualitative agreement with the analysis of DZ stars by 
Dufour et al.l (120071 ) who suggested that the presence of hydrogen in cool DZA white dwarfs 
could be explained better as the outcome of convectively mixed DA stars rather than due to 
accretion from the interstellar medium. This suggestion is based on estimates of the total 
mass of hydrogen present in the mixed hydrogen and helium envelope of DZ and DZA stars. 
Measurements of the hydrogen abundances in these stars suggest hydrogen l ayer masses in 
about the same range as that inferred here (see Fig. 13 of IDufour et al.ll2007l ). 



6. CONCLUSION 

We have presented a deta iled s t udy o f the spectral evolution of cool white dwarfs fol- 



lowing the pioneering effort of ISionl (I1984J ) . In particular, we have determined the ratio of 
helium- to hydrogen-atmosphere white dwarfs below T e g = 15, 000 K as a function of ef- 
fective temperature by increasing the size of the sample by a factor of 2, and by obtaining 
more precise measurements of the effective temperature of these stars through model atmo- 
sphere fits to Johnson V (or Stromgren y) and 2MASS JHKs photometric observations. 
We have confirmed that in this temperature range, the evolution of hydrogen- and helium- 
atmosphere white dwarfs is not independent, and we have revealed in greater detail the 
role of the convective mixing process responsible for the coupling between these two white 
dwarf types. We have found that about 15% of cool hydrogen-atmosphere DA white dwarfs 
between T e g = 15, 000 K and 10,000 K are transformed into helium-atmosphere non-DA 
white dwarfs at lower temperatures. Using a basic model for convective mixing, our results 
imply in turn that 15% of the DA stars that have survived the DA to DB transition have 



hydrogen layer masses in the range Mn/M tot ~ 10~ 10 to 10~ 8 , although the exact values 
depend on the assumed convective efficiency in the evolutionary model calculations. The re- 
maining DA stars would presumably have more massive hydrogen layers, probably in excess 
of M H /M tot ~ lO" 6 . 

Further analysis of the convective mixing problem would require a significantly larger 
photometric sample, such as the ugriz photometric sample of the Sloan Digital Sky Survey 
(SDSS). This sample covers a few areas of the sky with a completeness of around 95% 
up to g = 20. However, most spectroscopic fol l ow-ups have relied on c olor-color diagrams 
to identify white dwarfs (IKleinman et al.l 12004 lEisenstein et al.l l2006bl ) . and this method 
is strongly biased towards hot white dwarfs (T e g > 12000 K). Therefore, the SDSS white 
dwarfs identified so far cannot be used to study the convective mixing process at low effective 
temperatures. One way to pro ceed to identify co oler white dwarfs in the SDSS is through 
proper motion-color diagrams. lHarris et al.l (120061 ) have used this method to identify nearly 
6000 white dwarfs but a complete spectroscopic analysis similar to that of iKilic et al.l (120061 ) 
has yet to be performed to take full advantage of this sample. 
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Table 1. Comparison of Effective Temperature Determinations 







T cB (K) 


T ctt (K) 




WD 


Spectral Type 


This study 


Others 


Reference 


0046+051 


DZ 


6716 


6220 


2 


0135-052 


DA 


7353 


7275 


1 


0310-688 


DA 


15891 


15500 


1 


1142-685 


DQ 


8429 


7900 


3 


1647+591 


DAV 


12557 


12258 


1 


1917+077 


DBQA 


11488 


10200 


1 


1940+374 


DB 


15565 


16877 


5 


2032+248 


DA 


20287 


19511 


1 



Note. — R eferences: (1) A. G iannin as (2007, priv a te co mmu- 
nication), (2 ) iDufour et alj <2007h. (31 [Dufbur et al.l d2005h , (4) 
lOswalt et all lll99ll) , (51 iBeauchamp et al.l dl996fh 
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Table 2. Ratio of He- to H-atmosphere white dwarfs 



T eff range 



(10 3 K) 


^Hc-atm 


■^H-atm 


Ratio a 


15-13 


11 


41 


0.268 


13-11 


9 


38 


0.259 


11-10 


10 


44 


0.265 


10-9 


12 


41 


0.347 


9-8 


14 


38 


0.444 


8-7 


19 


49 


0.485 


7-6 


18 


51 


0.454 


6-5 


14 


38 


0.448 



a The ratio is corrected for the volume effect 
discussed in the text. 
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Fig. 1. — Location of the hydrogen convection zone (hatched region) as a function of effective 
temperature in the pure hydrogen envelope of a 0.6 M & DA white dwarf calculated with the 
ML2/a = 0.6 version of the mixing-length theory (from G. Fontaine & P. Brassard 2006, 
private communication). The depth is expressed as the fractional mass above the point of 
interest with respect to the total mass of the star. The thick solid line corresponds to the 
photosphere (tr ~ 1). 

Fig. 2. — The rati o of non-DA to DA white dwarfs as a function of absolute visual magnitude 



for the sample of ISionl (119841 ) . The error bars correspond to statistical uncertainties. The 
temperature scale for pure hydrogen model atmospheres with M = 0.6 M & is shown for 
comparison at the top of the figure. 

Fig. 3. — Predicted V magnitude as a function of T e g for pure hydrogen (solid line) and 
pure helium (dashed line) model atmospheres at log g = 8 and for which i^MASS = 16.5. 
The turnover of the hydrogen-atmosphere curve below 5000 K is due to the onset of the 
collision-induced absorptions by molecular hydrogen. 

Fig. 4. — Difference in (V — J) color indices as a function of T e g between logg = 8.5 
and 7.5 model atmospheres with pure hydrogen (solid line) and pure helium (dashed line) 
compositions. The general behavior of these curves can be explained in terms of the various 
opacity sources that peak at slightly different effective temperatures for these two logg values. 

Fig. 5. — Sample fits to the observed VJHKs photometric observations (error bars) of eight 
bright white dwarfs with hydrogen and helium atmospheres. The model fluxes are shown by 
open circles and the atmospheric parameters are given in each panel. 

Fig. 6. — Effective temperature variation, AT e g, for a corresponding change of the (V — J) 
color index of A(V — J) = 0.08 as a function of T eS for our pure hydrogen (solid line) and 
pure helium (dashed line) model atmospheres. The dotted horizontal lines represent the 
various bin widths of the histograms used in § 5. 

Fig. 7. — Differences in effective temperature determin ations based on our 2MASS sample 



(Tphot) and the photometric BVRIJHK data (T B rx) of iBergeron et al.1 (119971 . 120011 ) for the 
133 cool white dwarfs in common. The solid line represents a perfect match between both 
data sets, while the dashed lines correspond to temperature uncertainties equal to the bin 
widths of the histograms used in § 5. 

Fig. 8. — Top panel: Differences between our photometric T e g values and those determined 
spectroscopically for a common sample of 198 DA white dwarfs. The small open circles 
represent stars for which the 2MASS fluxes are close to the detection limit with formal 
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uncertainties in only two bands. The filled squares correspond to known double degenerate 
systems. The horizontal and dashed lines have the same meaning as in Figure [7J Bottom 
panel: Same as top panel but for 18 DB white dwarfs. 

Fig. 9. — Histogram of the number of hydrogen-atmosphere (single hatched) and helium- 
atmosphere (double hatched) white dwarfs from the 2MASS sample as a function of T e g. The 
number of stars in each bin is per unit of 10 3 K in T e g . 

Fig. 10. — The ratio of the observed space volume for helium-atmosphere and hydrogen- 
atmosphere white dwarfs in a if magnitude limited survey as a function of T e g. 

Fig. 11. — The ratio of the number of helium- to hydrogen- atmosphere white dwarfs as a 
function of effective temperature, corrected for the volume effect displayed in Figure \IU\ and 
discussed in the text. The filled circles use the binning defined in Figure [9] while the open 
circles present the results for the temperature bins shifted by 500 K. The error bars are 
statistical. 
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Figure [T] 
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Figure [3] 
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Figure H] 
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Figure [7] 
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Figure M 
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Figure [TO] 
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